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1. INTRODUCTION AND EXECUTIVE SUMMARY 

This report summarizes achie\ ements during year 3 of our project to in\estigate the use of ERS-1 
SAR data to study Arctic ice and ice/atmosphere processes. The project was granted a one year 
extension at no cost, and this report outlines our goals during the final year. The final report will be 
deli\ ered in October 1995. 

1 . 1 Project Objectix es 

The specific objects of the project are to: 

• determine how the development and ex olution of open xvater/thin ice areas xvithin the 
interior ice pack x ary under different atmospheric synoptic regimes; 

• compare how open xxater/thin ice fractions estimated from large-area divergence 
measurements differ from fractions determined by summing localized openings in the pack; 

• relate these questions of scale and process to methods of observation, modeling, and 
ax eraging over time and space; 

• determine xxhether SAR data might be used to calibrate ice concentration estimates from 
medium and loxx -rate bit sensors (AVHRR and DMSP-OLS) and the Special Sensor 
Microwax c/Imager (SSM/I); and 

• inx estigate methods to integrate SAR data for turbulent heat flux parametrization at the 
ocean/atmospherc interface xvith other satellite data. 

1.2 Year 3 Progress 

Major accomplishments ox er year 3 included: 

• The SAR-based ice classification xxork performed in the prex ious txvo years has been 
completed. The general effectiveness of SAR ice type classification xx'as assessed by 
companson xx ith Landsat TM data, and a specific evaluation of the performance of the ice 
type classifier used operationally by the Geophysical Processing System (GPS). 
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• Major insights have been gained into the seasonal cycle of typical pack ice and how that 
cycle is reflected in SAR imagety. B} a detailed qualitative examination, and the 
development and evaluation of a simple numerical/backscatter model, we are able to identify 
major features of the seasonal cycle and describe their appearance in SAR data. 

• As a consequence of the above, a time senes of mean SAR backscatter through the year 
was used as a means of evaluating both other remote sensing products (SSM/I) and a 
commonly used physical model of ice processes. 

• A simple methcxl for using SAR data for determining ice concentration was developed. 
By taking the effects of wet snow and meltponds into account, this method can be used to 
estimate the fractions of different ice types throughout the year. Although no in situ data 
vv as av ailable for v alidation, the concentrations obtained from this method appear reasonable 
when compared with those from other sources. 

• A study of the Ladv Ann Strait p»olynya using SAR imagery was conducted in which we 
learned how to interpret SAR imagery in the dynamic and sometimes confusing polynya 
environment and examined the change of backscatter with time on both short and long 
lempcvral scales. 

• In conjunction w ith a separate NASA grant, SAR imageiy and SAR-derived ice motions 
were combined with AVHRR-derived ice motions, SSM/1 data and sea-ice model output to 
siudv the relationships between observed ice concentrations and simulated divergence and 
shear dunng October 1991 and extended through May 1992 for the Beaufort Sea. No large 
dcscrepancies w ere found between simulated changes in concentration and those observed in 
the imagery. How ev er, the SSM/I data suggest substantially larger increases in first-year ice 
fraction, which might be related to sub-time-step changes in open water coverage. 
Significant coirelations exist between changes in SSM/I-derived ice concentration and 
simulated div ergence and shear over abevut 40% of the study area. 

1.3 Year 4 Plans 

During year 4 we intend to pursue the following activities: 

• The studv (M' the seasonal cycle of pack ice using SAR data will be continued. A particular 
objectiv e will be to vary the assumptions used in the numerical/backscatter model. Further 
efforts will be made to use the SAR time series as a tcxil for evaluating models and data from 
other sensors. 
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• The Lady Ann Strait and North Water polynya studies will be pursued further. Objectives 
will be to determine the typical evolutionary sequences of the polynya as a whole 
preparatory to a comparison with predicted ice growth and dynamics from a coupled 
ice/ocean polynya model. 

• Comparison of SAR and AVHRR-derived ice motions, SAR lead patterns, SSM/I ice 
concentrations, and model simulations will be expanded to include the period from June 
1992 - June 1993, in conjunction with separate NASA funding. 
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2. YEARS ACTIVITIES 


2. 1 Ice type classification 

Dunng the past year, we improved and consolidated our previous work on SAR-based ice type 
classification using Landsat TM data as "truth". A more careful seleetion of training data regions 
was made and the coregistration, calibration and statistics prtx:essing improved and streamlined. 
Our stud\ of ice typing, which has been described in a paper currenth in press (Steffen and 
Heinrichs 1994a), had two major objectives. The first was to evaluate SAR-based ice typing in 
general using Landsat TM data, and the second and more specific was to assess the performance 
of the operational classifier used in the Geophysical Processing System (GPS) at the Alaska SAR 
Facility (ASF). 

The data set used in our work consisted of two o\ erlapping pairs of Landsat TM and ERS-1 SAR 
images taken on April 15 and 18, 1992 covering a portion of the Beaufort Sea centered at 145® W, 
73® N. The Landsat data was obtained from EOSAT Corporation and was geolocated, mapped to 
the UTM projection, and gridded to 25 m spatial resolution (slightl) better than the aetual 30 m 
resolution) in order to better mateh the SAR. Data from the C-band (5.3 GHz) SAR instrument 
aboard ERS-1 was obtained from the ASF at its full 12.5 m spatial resolution, resampled to 25 m 
resolution, and mapped to the UTM projection. Coregistration of the image pairs was a 
challenging process, gi\ en both the movement of the ice and certain errors in the data processing 
ol the SAR, and w as accomplished by the use of successively smaller offsets determined visually 
and the use of color composite images for precise registration. 

Once coregistration and image subsetting were complete, images of TM band 4 (760-900 nm) 
and, (x-casionally, TM band 6 (10.4 - 12.5 mm in the thermal infrared) were used to select 
training data regions lor each of seven ice types as defined by the standard WMO nomenclature. 
We accumulated training data areas until a sufficient sample size (roughly 1000 points) w'as 
attained for each type. The training data were used to define a super\ ised classification (maximum 
likelihotxl) based on Landsat TM band 4. 

Table 1 shows the results of our statistical analysis of the SAR pixels assigned to eaeh ice type as 
defined b} the Landsat supervised classifieation. A graphical representation of these results is 
shown in Figure 1. The TM data shows a distinct separation betw'een most ice classes w'ith, 
howc\er, some o\ erlap between first year ice (FYI) and old ice (OI). This is due to the snow' 
co\er over both ice types causing their albedos to be very similar. The SAR data, however, 
show s major ov erlaps between many ice types. Ol and FYI can be readily distinguished in the 
SAR data, but the nilas (NI), grey ice, and grey-white ice categories have substantial overlap, 
both with each other and with the OI and FYI classes. The ice free (IF) areas also have 
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backscatter \ alues similar to the solid ice types, but dilTerent lor the iwo days studied. This i*s a 
consequence of wind roughening (increased backscatter caused by wind-formed capillary waves 
on the cx:ean surface). Wismann (1993) has provided a fomiula which relates the backscatter of 

IF areas backscatter to wind speed (cK^n = w here is the backscatter of the 

ice tree surface and u is the wind speed). The c'lbseiwed backscatter from ice free areas is slightly 
less than w hat w ould be expected from this relationship, but dc:>es show; a higher backscatter on 
April 18 when winds were stronger. We hypothesize that the Wismann formula is more 
appropnate for larger open ocean areas then the leads found in the pack ice environment. 

[able 1. Derived SAR backscatter coefficient (ct°) statistics (mean mid standard de\ iation) by ice type 
based on the supervised ice t>pe classification from I^idsat Thematic Mapper chamiel 4 for the 
Beaufort Sea (April 16 and 18, 1992). 



April 16, 1992 

Backscattermg Coefficient (dB) 

April 18, 1992 

Backscattering Coefficient (dB) 

Ice Type 

Mean 

SD 

Mean 

SD 

IF 

-15.12 

3.28 

-5.31 

1.95 

NI 

-12.35 

4.03 

-10.39 

3.95 

GI 

-14.01 

2.98 

-13.14 

2.38 

GWI 

-15.17 

3.53 

-14.25 

2.99 

FYI-I 

-16.26 

4.67 

-16.75 

3.23 

FYI-n 

-17.41 

4.87 

-16.70 

3.69 

OI 

-10.97 

3.43 

-10.20 

2.47 


i igurc 1 Graph showing meiuis and standard de\ iations of Gindsat TM band 4 six^ctral radiance and 
s AR backscatter coelficient lor each of seven ice types. 


ICE TYPING: Apri 18. 1992 

LANDSAT THEMATIC MAPPER CHANNEL 4 (760-900 NM) 
ERS-1 SAR (Otend) 


ICE TYPING: April 16. 1992 

LANDSAT THEMATIC MAPPER CHANNEL 4 (760-900 NM) 
ERS-1 SAR (C4Mnd) 






Our evaluation ol the GPS classifier was based on the results of the SAR/Landsat comparison 
described above. The GPS ice types are somewhat different from the WMO list (Kwok and 
Cunningham 1992, Kwok et al 1992), and the training data was recombined to reflect this 
difference. Table 2 lists the SAR means and standard deviations for each of the categories as 
defined by the GPS and the TM supervised classification for April 16. While, the two 
classifications agree well lor PWl-1 and Ol, significant differences in the means exist for the 
combined IF/Nl/Yl and the FYl-Il classes. The GPS threshold used for IF generally restricts 
identification to calm wind conditions. IF areas identified in the TM imageiy, because of the 
moderate winds on April 16, have higher backscatter v alues and are misclassified by the GPS as 
other types. Bv' calculating the percentage of correct decisions made by the GPS classifier as 
compared to the TM-derived classification, we were able to quantitatively assess the GPS 
perfoiTnance. These results are listed in Table 3. As expected from the means and standard 
dev iations, Ol classification accuracy is the best of the four GPS types with FYI-1 a distant 
second. IF/Nl/Yl and FYl-11 categories are clearly not performing well, most likely due to 
misclassification of IF pixels. Ice concentrations were determined for both study dates using SAR 
and TM thresholds. The TM values are closer to the expected ones, and show the expected 
relation to divergence (lower concentration on April 18 asscxiated with higher wind speeds). 
SAR-derived concentrations show the opposite tendency, most likely due to a larger fraction of IF 
pixels being wind-roughened and misclassified as solid ice. 

fable 2. .Mean S AR backsaitter coefficient iuid standard deviation (SD) for the GPS ice classifier and 
the l.andsal TM chiumel 4 based ice classification for ice free, nilas and voting ice (IF, NI, YI), thin 
lirsl-vear or smooth first-yeiu" ice (FA'I-I), metliuin and thick first-year or deformed first-year ice (F\'I-II) 

:uid old ice (Ol) on .April 16, 1992. 


Ice Type 

GPS mean 

GPS SD 

TMmean 

TMSD 

IF, NI, YI 

-20.7 

3.5 

-14.2 

3.5 

FYI-I 

-17.8 

3.9 

-16.3 

4.7 

FYI-n 

-14.0 

3.7 

-17.4 

4.9 

Ol 

-10.8 

3.3 

-11.0 

3.4 
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Table 3. Perl'onnance of the GPS SAR ice tjpe classifier. The correctly classified percentage for each 
SAR ice type as compared to tlie TM ice classes is given in parentheses. Tlie comparison was done 
with the co-registered I^idsat smd SAR images (650,000 pixels) for April 16, 1992. The total 
percentage of each ice type as derived from TM for the entire image is also given. 


SAR 

Landsat TM Channel 4 Ice Qasses 

Ice Classes 

IF&YI 

FYI-I 

FYI-n 

OI 


(0.3 %) 

(20.8 %) 

(21.1 %) 

(57.8 %) 

IF, NI, YI (%) 

(13.6) 

34.9 

31.2 

20.3 

FYI-I (%) 

11.9 

(36.1) 

27.1 

24.9 

FYI-n (%) 

40.5 

21.0 

(16.4) 

22.1 

OI (%) 

4.4 

7.2 

25.8 

(62.6) 


The final portion of the study in\olved an examination of the effectiveness of using the SAR and 
TM data together in a multispectral classification. SAR is superior for separating FTI and OI, 
where Landsal has advantages for IF, NI, and YI. We assembled a set of training data using 
color composite images of SAR and TM data together to assign a priori categories and evaluated 
the performance cM se\ cral channel combinations against the training data. For each channel set 
in\ e.sti gated, we calculated the parameters for a maximum likelihood classifier, which was then 
applied to the training data. Accuracy percentages for each category’ were calculated for each 
channel set. The combination of SAR and single TM visible channels (channels 2,3, or 4) was 
lound to improve accuracy over TM for all but the FYI types. Adding multiple visible channels 
impro\ed accuracy for FYI-II, and adding channel 6 (thermal infrared) resulted in an accuracy 
greater than or equal to 90 % for all types. 

2.2. Studies of the Seasonal Cycle of Sea Ice Using SAR Data 

The potential of SAR data for longer-term monitoring of sea ice has only begun to be explored, 
despite the important ad\antages of SAR data for this kind of research. We conducted a study 
aimed at determining the \ ariation of backscatter with time from typical pack ice and how these 
changes relate to physical changes in the ice are important problems 

The data used in this study covered a region in the Beaufort Sea centered at 78° N, 153° W during 
a 13-month period from October 1992 to November 1993. A series of 31 images from the ERS-1 
C-band SAR was obtained from the Alaska SAR Facility (ASF) and calibrated using the 


Page 8 




operational ASF algorithm. As part of the calibration process, the data was scaled from Digital 
Number (DN) values to backscatter \alues in dB. In addition to the SAR data, we obtained 
concurrent meteorological data from the National Meteorological Center (NMC), including surface 
air temperature, pressure, and geostrophic wind fields. The complete time series is illustrated in 
Figure 2. 


f igure 2. Graph ol the inemi ERS- 1 S.\R backscatter <uk1 N.MC temperature juicl wind speed data used 
to study the seasonal cycle of sea ice in the Beaufort Sea. Ronuui numerals identity the regime (see 
te.xt for e.xplanation) and \ ertical lines are regime boundaries. 



1991 


1992 


2.2. 1 Qualitative Examination 

The first general obsen ation we made is that there is substantial \ariability in the appearance of 
sea ICC throughout the year in SAR imagery. Not only may a given value of backscatter 
conespond to more than one kind of ice (or e\en wind-roughened IF), but this relationship 
changes through the year. Nevertheless, we have segmented the year into three major 
c\ tilutionary regimes that, we feel, have a consistent set of ice signatures and a correspondence to 
major features in the seasonal cycle of sea ice. 

The first regime w c identified (regime 1) cov ers the majority of the year, beginning in autumn and 
extending all the wav through spring until summer melt occurs. The signatures of FYI and 
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multiyear ice (MYI) are relatively stable, about -16 and -10 dB, respectively. The major 
modulation of the mean backscatter in the time senes was the \ ariations in ice concentration and 
the w ind roughening and corresponding backscatter increases of IF. Beginning shortly before the 
summer solstice we have another major regime (regime II), in w'hich the appearance of sea ice 
changes dramatically. First, the snow cover above the ice surface begins to melt. At first, this 
in\ olves primariK the addition of w'ater to the snow layer. This water, because it is embedded in 
the snow matri.x, is not subject to wind roughening and acts as a v ery efficient absorber of 
electromagnetic radiation at radio wavelengths. Because of this, the MYI and FYI have a 
significantly reduced backscatter. Pressure ridges, w'hich protrude above the layer of wet snow, 
appear as bright linear features. Later in the summer, the wet snow melts completely, w'ith much 
of the water collected in meltponds. The meltponds are also subject to w ind roughening, and the 
result is a dramatic increase in the sensitivity of the mean backscatter to wind speed. The floes 
also become disaggregated. At the end of the summer, the appearance of sea ice in the SAR 
imagery changes once again, defining the third regime (regime III). MYI regains its normal 
appearance as the meltponds refreeze, and FYI begins to form once again. The major difference 
in this regime from the pre\ ious regime is the presence of broadly linear bright features. These 
features have a general appearance similar to that of the IF areas in the summer, and we suspect 
the\ are due to the refreezJng of the leads with their multitudes of disaggregated floes combining 
in a jumble of angular fragments w hich act as reflectors (both direct and ^'ia cornering) for 
radiation at C-band wavelengths. 

2.2.2 A Simple Ice Co\ er/backscatter Model 

In order test our hyp(')theses and the conceptual model that originated during the qualitative 
examination of the data set, we constructed a simple numerical model. Our plan was to link a 
simple backscatter model to this numerical model and then compare the resulting backscatter 
estimates to the actual time series. 

Wc began by assuming ice concentrations to be constant in each regime. This assumption 
simplifies the initial e\aluation of the model and, furthermore, can be changed at any time. The 
physical charactenstics of the model are, as would be expected, dependent on the time of year 
(using the regimes defined above). In regime I (autumn/winter/spring), we have only FYI, MYI, 
and IF present (NI, grey ice, and grey-white have been excluded to simplify the model). During 
the summer, all of the FYI becomes MYI, which is partially co\ ered by wet snow and meltponds. 
Because of their importance to the energy balance of the ice, and because of their strong effect on 
the SAR backscatter, we added a more detailed parameterization of their growth. A maximum 
fraction for snow co\ er on the ice is chosen, as is a maximum meltpond fraction. These numbers 
were assumed to be 0.65 and 0.4, respectively. At the beginning of the summer season (defined 
b\ the excursions of the surface air temperature above and below 0° C) the entire maximum 
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1 faction ol wet snow appears. During the summer, the meltpond fraction is calculated for each 
image by a simple linear relationship of the day since the beginning of summer and the maximum 
meltpond fraction. Since the meltpond formation is at the expense of the snow -covered area, the 
area of snow on a gi\en day is reduced to the maximum meltpond fraction linearly with the 
number of days since the beginning of the melt season. During the late summer freeze (regime 
III), no IF is allowed and all of the IF from regime II is assigned to the bright ice (BI) category. 
Table 4 lists the ice concentration values used for the different regimes in our ice co\ er model. 

Table 4. Surface type fractions (range 0.0 to 1.0) used in ice cover model by regime 


Surface l ype 

Surface Type Fraction 

Regime 1 

Regime II 

Regime III 

MYl 

0.85 

0.9 

0.85 

fYi 

0.1 

0 

0.1 

Bright Ice 

0 

0 

0.05 

If 

0.05 

0.1 

0 

Wet Snow 

0 

varies 

0 

Meltponds 

0 

varies 

0 


On lop ol this simple model of ice cover through the year we added an equally simple model of the 
mean backscatter of the .scene. The mean backscatter of the scene is, by definition, the average of 
the backscatter \alues at each pixel. By segmenting the pixels into types, we can calculate the 
mean backscatter by adding the backscatter of each type weighted by the fraction of area covered 
by that type (i.e., the concentration). Including all of the ice surface types described above, we 
then ha\'c; 


+Aifo0if + AbicV^’bI + AMpOt»Mp + AwsO^W’S (1) 

where ^ <^Bb ^nd a0\vs i^re the backscatter values I'or multiyear ice, first- 

year ice, ice free, bnght ice, meltponds, and wet snow respectively. Amy, Ar-, Aif, Abi, Amp, 
and A\\ s are the corresponding cover fractions for each lypx;. We assumed that each of the solid 

icc types has a characteristic backscatter, and that all pixels of that type have that backscatter value. 
To obtain these \alues, rather than the rigorous training data selection procedure used for the ice 
type classification study, we examined the time series and chose what we considered typical 
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values for each type = - 10.5 dB, o0^T^ - - 16.0 dB, o<^>b[ - -7.0 dB, and = -22.0 dB). 
We assumed that open leads and meltponds are both subject to wind roughening, so that and 

o^’iF are identical and can be calculated by the Wismann formula described above using the NMC 
w ind speeds. It was assumed that the wind speed is constant ov er the scene. 

Figure 3 show s the results of running the combined ice cov er/backscatter mcxiel. It can be readily 
seen that the model reproduces the major features of the seasonal cycle, including the stable 
backscatter values in the winter and spring, the sudden backscatter drop in summer, and the rise 
during freezeup. More quantitatively, the RMS difference between the actual and simulated 
backscatter time series is 1.01 dB, which represents about 6 % of the full range of backscatter 
values obsen ed in nature.Some differences exist, however, particularly on shorter time scales. 
The simulated backscatter seems to be out of phase with the actual time series during the spring 
months. This may be a function of short-term variations in ice concentration due to large-scale ice 
motions or the passage of spring storms and consequent divergence. 

Figure 3. Comp^irison of actual SAR memi backscatter time series with simulated time series from 
simple ice co\er backscatter model, SSM I NAwSA team algorithm, and Hibler ice 
d>Tiamics thermod}iiamics model. 


ACTUAL vs. SIMULATED BEAUFORT SEA SAR TIME SERIES 



Oct Nov Dec Jan Feb Mar Apr May Jun Jul Aug Sep Oct Nov Dec 
1991 1992 

Month 

ACTUAL MEAN SAR BACKSCATTER 

SAR BACKSCATTER SIMULATED WITH SIMPLE MODEL 

SAR BACKSCATTER SIMULATED WITH SSM/I ICE CONCENTRATIONS 

SAR BACKSCATTER SIMULATED WITH HIBLER ICE CONCENTRATIONS 
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In order to lest our simple model, as well as to in\ estigate the use of the backscatter time series as 
a \ alidation tool in general, we obtained ice concentration estimates for the study area from two 
other sources. The first was the SSM/1 sensor, using the NASA team algorithm with global tie 
points. The second source was a Hibler ice dynamics/thermodynamics model with viscous-plastic 
rheology (Hibler 1985), run for the entire Arctic, and from which we extracted the data for the 
gnd point nearest the study area. The Wismann relationship for the backscatter of IF and the 
backscatter values lor the solid ice types were used to calculate mean backscatter from each ice 
concentration time series. Being a unique feature of our ice cover model, we did not include the 
effect of meltponds and wet snow in the summer. 

Figure 3 also show s the actual SAR time series compared with the simulated SAR data obtained 
from all three sources (our simple model, the SSM/1 NASA team algorithm, and the Hibler 
model). The SSM/l-denved backscatter estimates are clearly too low, reflecting the generally 
underestimated concentrations that have been obserxed using global tie points (Steffen and 
Schweiger 1990, 1991), but tend to track the changes in the real time series fairly w'ell. The 
Hibler-denx ed estimates, on the other hand, have values of about the same overall magnitude as 
the actual SAR measurements in the winter and spring, but show Uxi little \ ariation all year and 
particularly in the summer. Table 5 shows the results of a statistical comparison between all three 
simulated SAR time series and the actual measurements. The lowest RMS error is from the 
SSM/1, which as mentioned tracks the real time series fairly well, particularly in summer. The 
Hibler and our simple ice co\ er models haxe about the same RMS error (slightly higher than that 
ft>r the SSM/1), surprising considering the relative complexity of the two. The mean 
actual/predicted difference is greatest in magnitude for the SSM/I, confirming the observation of a 
significant bias in the estimates. The mean absolute difference, perhaps the most robust measure 
ol statistical difference in cases where the distribution if non-normal, is almost identical for all 
three. 


l able 5, r.rror statistics I'or actual \s. simulated mean SAR backscatter (in dB) for tlaee different 
sources of simulated SAR. 


Source 

RMS Error 

A\ erage l irror 

Average Absolute 
Error 

Cl ice cover backscatter 
model 

1.01 

-0.30 

0.69 

SSM 1 N AS/\ temii 
cdgoritlim 

0.69 

0.51 

0.66 

Hibler ice model 

0.96 

-0.22 

0.64 
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2.2.3 Ice concentration Algorithm Development 


Estimates of ice concentration are useful for both climate research and operational applications 
(shipping, in particular). During the summc months, when melt cxcurs and the ice surface 
undergoes its most dramatic changes, estimates from passive microwave sensors are notoriously 
inaccurate (to the extent that they are not calculated operationally). The time series analysis 
described above suggested the a possible approach to obtaining ice concentrations throughout the 
year from SAR data, with the availability of the NMC meteorological data to help constrain the 
problem. 

At an\ time the fractional coxerage of the \ arious surface types must add to unity. Furthermore, 
the sum ol the backscatter from all sources vxeighted by the concentration of that ice type must 
equal the mean backscatter of the scene. Thus we have, assuming the presence of only FYI, 
MY I, and IF, we have 


Amy + Apv +Aif =1 (2) 

cTfO = Am\xY\p^' + Ap^-cY^FY +Aif (3). 

B\ soh ing I'or Ajf and substituting (2) into (3), we can obtain an equation relating Ap\’ and Am\- 
From our prev ious work, and o^^fy ttrc known (cY^ix = 0 in summer). From the NMC wind 

speeds, and the Wismann formula, o*^if can be calculated. 


Wc begin calculation (vl‘ icc concentration by counting the number of pixels with backscatter values 
above and below a FYI/MYI threshold (-13.5 dB) and obtain AFX(measured) and AMY(measured)- 
Because the entire range of pxissible backscatter v alues has been divided into FYI and MYI, the IF 
points must fall into one of the two classes. Since cY*]i is known, wc can decide in which ice 

type the water points have been included. The other type can then be assumed to be free of IF 
points and thus exact (i.e., Apv- = AFY (mea.sured) or Axn = A.M^Ymeasured))- Since we have an 
equation relating Afy and A\p , we now know both, and calculation of A]f is a simple matter of 
substitution. 

Figure 4 shows the icc concentrations calculated using this method compared with those from the 
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SSM/I sensor and Hibler model. It can be obser\ed readily that SSM/I concentration values for 
MYI are too low , and dramatically so in the summer, when operational ice concentration estimates 
are not produced. In general, the SAR-deri\ed concentrations are closer to those obtained from 
the Hibler model than the SSM/I, and certainly fall in the reasonable range. A detailed statistical 
exaluation has not been performed so far, due to the preliminaiy and exploratory nature of this 
work. Potential problems with our approach include eixors w hen the wind speed is close to the 
threshold x alue used to separate FYI from MYI (perhaps not significant gix en the veiy narrow 
range of wind speeds yielding such backscatter \ alues). 

Figure 4. Ice eoncenlralioiis obtained from S.\R algorithm, N’,\S.\ team algorithm (SSM 1) and Hibler 

model. 



Multiyear Ice 

■ First Year Ice 

Ice-free (Open Water) 
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2.3. Lady Ann Strait Polynya Study 


It is becoming more widely realized that coastal polynyas ha\c a major inlluence on the heat and 
energy balance of the ice-covered seas. Furthermore, the large rates of ice fonnation and brine 
prtxluction typical of polynyas are important climatic factors. Our objecti\ e was to e.xamine ice 
co\er \ ariations and the underlying physical processes in one polynya, making use of the rapid 
repeat coverage during the ERS-1 ice cycle and the SAR instrument's high resolution and ability 
to penetrate clouds and darkness. Because of extensive previous study, and high suitability for 
SAR observation, we chose the recurring polynya located in Lady Ann Strait at the mouth of 
Jones Sound in northwestern Baffin Bay as our subject. 

We obtained, from ASF, a series of 30 images of the Lady Ann Strait region covering a period of 
5 months between Now 20, 1991 and April 17, 1992. The images were of the low resolution 
prcxluct, with 100 m resolution, each 100 km on a side. The images were calibrated using the 
standard GPS technique, converted to decibels (dB), and scaled. 

As can be seen from Figure 5, the appearance of the polynya goes through dramatic changes on 
short time scales. The Feb. 6 image shows the polynya filled with brash ice composed of broken 
new icc. Three da\ s later enough ice has formed to co\ er the polynya, and by the next image this 
ccn er has broken up as w ell. A few relatively large YI floes are \ isible. By Februaiy 15, all of 
this material has been removed by wind action and ice is again beginning to form. A more 
detailed understanding ol the kinds of ice visible in the SAR images was obtained by taking 
transects through selected areas of the polynya and examining the backscatter values along the 
transects. One such profile (taken from the February 6 image) is shown in Figure 6. From this 
profile, as w ell as others, w e have noted an approximately linear increase of backscatter with ice 
thickness for new ice and NI. 
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Figure 6. Backscatter profile taken across portion of l^dy Ann Strait Polynya on Febniarv’ 6, 1992. 




B\ taking profiles across stable FYI and MYI through the duration of the study period, we were 
able to examine the temporal change of backscatter with time over longer time scales. Figure 7 
show s the results of this analysis. The MYI floes embedded in the FYI matrix show virtually no 
change in the backscatter, and neither does the smooth FYI. For FYI in general, however, as well 
as new ice at the ice arch, a distinct downward trend in the backscatter is observ ed. At present, we 
do not ha\ c a ngorous explanation for these trends. 


I igurc ^ Sea ice backscatter time series derived from HRS-1 SAR profiles for the lady .Aim Strait 
)X)lynya. 



1992 
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The first conclusion that we can draw from our work on the Lady Ann Strait polynya is a negative 
one. We do not feel it is possible to obtain ice motions within polynyas from ERS-1 SAR data. 
The ice orbit of 3 da>s does not allow ice dynamics to be resohed. Ice velocities within the Lady 
Ann Strait polynya frequenth reach 20 km/day, and can attain 50 km/day (Ito, 1981). In a 3 day 
period, the ice within the polynya can travel 60 -150 km and move' completely outside the 
bcvundaries of a 100 km square ERS-1 scene. Radarsat should prov ide sufficient temporal 
resolution to allow accurate measurement of polynya dynamics. 

Despite the difficultv in resolving ice motions, the SAR data has allowed us to observe the 
variability of ice cover in the polynya over a relatively long time pericxl. Although much 
inteipretation of the SAR imageiy is required in order to assess the ice conditions within the 
polynya, the high spatial and temporal resolution of our data set makes it a v aluable resource. By 
examining backscatter values along transects within the polynya, we have gained further 
understanding of the way in which backscatter v aries with ice growth for veiy young ice. 

The work described above is a vital step towards the ultimate goal of being able to monitor the ice 
cov er within polynyas and thus determining their heat and energy balances and influences at local, 
regional, and global scales. Techniques have been developed (ice classification by profile, for 
example) that w ill be invaluable for polynya study when more suitable data sets become available. 
A paper describing our study of the Lady Ann Strait polynya is in preparation and will be 
submitted .shortly to the Journal of Geophysical Research (Steffen and Heinrichs 1994b). 

2.4 Response of the Ice Pack to Atmospheric Synoptic Systems 

Ice concentrations retneved from SAR, SSM/I, and AVHRR imagerv' and simulated using a 
dynamic-thermodynamic ice model were compared in relation to the passages of strong 
atmospheric low-pressure systems in the western Arctic for a case study in October 1991, and for 
a lime series through May 1992. Combining different image types and ice model output to study 
ice conditions allows us to minimize some of the uncertainly in the different data types, and carries 
more confidence in the interpretation of actual conditions. Since the term "ice concentration" can 
represent different conditions as defined for different sensors and models, care is needed when 
intercomparing data sets. 

With this caveat in mind, magnitudes of changes in ice concentration as depicted in SSM/I and 
model output agree reasonabl}’ well in the case study for the October (including AVHRR) and 
over the autumn, winter, and spring period into summer 1992 (SAR, SSM/I, and model output 
only), showing a 3-5% reduction in ice concentration during the passage of strong storm systems. 
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SAR imagen shows some e^ idence of increased lead de\ clopment and/or opening during passage 
of the systems, but the total increase in IF/thin ice fraction amounts to onh an additional 1 to 2%. 
AVHRR imageiy show greater numbers of leads than are seen in the SAR data, but none of the 
data sets indicate a large-scale disruption of the ice co\er such as has been seen in late summer 
cases (e.g., Barry and Maslanik, 1989). The largest changes in the ice pack are indicated as an 
increase in first-year ice fraction in the SSM/1 data and model output, which is consistent with 
refreezmg of newly-formed IF areas under di\erging and/or shearing conditions. All the data 
types agree within a few percent for means of total concentration during the autumn through 
spring penod. Disagreement is substantial among the estimated proportions of first-year and 
multiyear ice. 

Differences in simulated concentrations between the runs using viscous-plastic and cavitating fluid 
rheologies results from the apparent tendency for the ca\ itating fluid rheology (used here without 
a representation of shear strength) to o\erestimate drift speeds relative to a t iscous-plastic 
rheology dunng periods of rapid ice drift. When averaged over all observations, the mean 
difference between the two rheologies is small. The drift speeds simulated using the viscous- 
plastic rheology arc closest to the buoy-measured and remotely-sensed drift speeds. 

Di\ergence rates estimated using optimalh -interpolated AVHRR and SAR-deri^ed motion fields 
arc similar and show considerable local \anabilit\ v\ithin a mean converging field. Simulated 
di\ ergencc rates oxerestimate convergence relatixc to the obserxations, although both obsened 
and modeled rates arc small. Lcx;alized changes in SSM/I-derived concentrations are typically 
consistent w ith regions of conx ergence and dix ergence in the pack. Significant correlations were 
found in some areas betxx een changes in SSM/I-derix ed concentrations and simulated divergence 
and shear ox cr about 40 % of the study area. 

Enough uncertainties exist regarding data sets, algorithms, and definitions that the absolute 
amount xif open-xxatcr prtxluction in relation to atmospheric synoptic systems cannot be 
determined unambiguouslx . The effects of such differences in retrieved concentrations, when 
considered in tenns of turbulent llux estimates, shoxx substantial differences depending on the 
data type and classification scheme used to retriex e ice concentrations. 

To further define ice-pack responses to atmospheric conditions, xvork is underway to extend these 
case studies in combination xvith surface and aircraft obserxations, and to build a long-term time 
scries of remotely-sensed ice concentrations, ice motion, and simulations using different modeling 
options. In addition to specifying the requirements of remote sensing systems and field programs 
lor sea-ice applications, such xvork should assist in the eventual assimilation of such merged data 
sets into operational ice models, and should aid in testing simplifications desirable for climate 
modeling applications. 
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3. YEAR 4 PLANS 


3. 1 Ice Type Classification 

We do not anticipate significant further work on comparing SAR and TM-derived ice type 
classification w ith SAR. Possible investigations, if time allows, include the use of texture and 
other spatial information, along with single and multiple channel images, in the ice typing process. 

3.2 Seasonal Cycle Study 

The seasonal cycle study is a major focus of our work, and will be continued. The assumptions 
used in the simple ice cover/backscatter model will be reexamined, and runs will be made w'ith 
different \ alues for some parameters (for example, the backscatter and tractional coverage of wet 
snow and meltponds). Results of our other work on the backscatter of new ice and YI will be 
included to study the effects of these types on the mean backscatter from the SAR scenes. We 
intend to investigate the effects of using image-dependent thresholds for ice types rather than the 
global thresholds currently in use. Finally, w e will explore the inclusion of a more complex sea 
icc backscatter model instead of the extremely simple backscatter model used so far. 

3.3 Ice Concentration 

The ICC concentration algorithm we have de\ eloped would benefit from further investigation. We 
plan to study the inclusion ol' the thinner ice types (Nl, grey ice, and grey- white ice) in the 
algorithm as well as the use of multitemporal information to further constrain the problem. The 
use of an estimated lead formation rate may be a promising approach to improvement. 

3.4 Polynya Studies 

We afso expect to continue our analysis of the Lady Ann Strait polynya data set. The cause of the 
ob.sen ed backscatter trends on monthly time scales remains an open question, and we will 
examine and evaluate possible mechanisms. The SAR "snapshots" of the polynya, while not 
suitable for a rigorous study of ice dynamics within the polynya, may be useful for validating a 
coupled ice/ocean polynya model currently under development by J. Heinrichs. 
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6. ACRONYM LIST 


ASF 

Alaska SAR Facility 

BI 

Bright Ice 

CU 

University of Colorado 

EOSAT 

Earth Observation Satellite Corporation 

ERS-1 

Earth Resources Satellite - 1 

FYl 

First-year Ice 

GPS 

Geophysical Processor System 

IC 

Ice Concentration 

IF 

Ice Free 

IDL 

Interactive Data Language 

MYI 

Multiyear Ice 

N1 

Nilas 

NMC 

National Meteorological Center 

01 

Old Ice 

SAR 

Synthetic Aperture Radar 

SIPS 

Spectral Image Processing System 

SSM/I 

Special Sensor Microwave Imager 

TM 

Thematic Mapper 

UTM 

Uni\ crsal Transverse Mercator 

WMO 

World Meteorological Organization 

Y1 

Young Ice 
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